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Abstract: We propose an inflationary scenario, M-flation, in which inflation is 
driven by three N x N hermitian matrices i = 1, 2, 3. The inflation potential of 
our model, which is strongly motivated from string theory, is constructed from $j and 
their commutators. We show that one can consistently restrict the classical dynamics 
to a sector in which the $i are proportional to the N x N irreducible representations 
of SU{2). In this sector our model effectively behaves as an N-flation model with 
3iV^ number of fields and the effective inflaton field has a super-Planckian field 
value. Furthermore, the fine-tunings associated with unnaturally small couplings in 
the chaotic type inflationary scenarios are removed. Due to the matrix nature of 
the inflaton fields there are 3A^^ — 1 extra scalar fields in the dynamics. These have 
the observational effects such as production of iso-curvature perturbations on cosmic 
microwave background. Moreover, the existence of these extra scalars provides us 
with a natural preheating mechanism and exit from inflation. As the effective inflaton 
field can traverse super-Planckian distances in the field space, the model is capable 
of producing a considerable amount of gravity waves that can be probed by future 
CMB polarization experiments such as PLANCK, QUIET and CMBPOL. 

Keywords: Inflation, curvature and iso-curvature perturbations, preheating. 



* amj ad@umich . edu 

t firouz@physics . mcgill. ca 

^ j abbari@theory. ipm . ac . ir 



Contents 



1. Introduction 2 

2. M-fiation scenario, the setup 3 

2.1 Truncation to the SU{2) sector 4 

2.2 Consistency of the truncation to the SU{2) sector 5 

3. Various inflation models resulting from M-flation sceucU-io 7 

3.1 Chaotic inflation 7 

3.2 Symmetry breaking inflation 8 

3.3 Inflection point inflation 9 

4. Mciss spectrum of modes in M-flation 10 

5. Power spectra in the presence of ^'j modes 12 

5.1 Linear perturbations 12 

5.2 Power spectra from Hubble-crossing to end of inflation 14 

5.3 Iso-curvature vs. entropy perturbations 15 

5.4 Curvature and iso-curvature perturbations power spectra for specific 
examples 17 

5.4.1 Chaotic infiation 17 

5.4.2 Chaotic infiation 19 

5.4.3 Symmetry breaking infiation 20 

5.4.4 Infiection point infiation 21 

6. End of inflation and preheating 21 

6.1 Particle creation during slow-roll infiation is not harmful to M-fiation 22 

6.1.1 Quantum production of ^r,im modes during infiation 22 

6.1.2 Back reaction on the infiationary dynamics 25 

6.2 Particle creation and the preheat scenario 26 

7. Motivation from String Theory 27 

8. Discussion 29 
A. Symmetry breaking inflation 31 



- 1 - 



1. Introduction 



Recent observations, specially the five years Wilkinson Microwave Anisotropy Probe 
(WMAP5) data [1] strongly support inflation as the theory of early Universe and 
structure formation. In their simplest forms, models of inflation are constructed 
from a scalar fleld, the inflaton fleld, which is minimally coupled to gravity. The 
potential is flat enough, so that a period of slow-roll inflation is achieved. These 
simple models of inflation, not only solve the problems associated with the stan- 
dard big bang cosmology such as the horizon problem, the flatness problem and the 
monopole problem, but also generate quantum perturbations which seed the struc- 
tures in the Universe. The slow-roll models of inflation produce perturbations on 
cosmic microwave background (CMB) which are almost scale-invariant, adiabatic 
and Gaussian which are in very good agreement with the current data. Theoreti- 
cally, on the other hand, inflation still remains a paradigm and one can construct 
many inflationary models compatible with the current data. There have been many 
attempts to embed inflation in string theory, for a review see [2]. 

Theories of multi-field infiation, in which one deals with more than one scalar 
field, have also been studied [3], for a review see [4]. In the multiple field infiationary 
models one can perform a rotation in the field space of scalar fields where the infiaton 
field is evolving along a trajectory while the remaining fields are orthogonal to it. 
These extra fields, like the inflaton itself, have quantum fluctuations which once 
stretched to super-Hubble scales can become classical and can therefore contribute 
to the power spectrum of iso-curvature as well as curvature perturbations, the details 
of which depends on the post inflationary dynamics and the reheating scenario. A 
speciflc possibility in the multi-fleld inflation is the idea of assisted inflation [5]: while 
the potential is too steep for an individual field to support infiation, the collective 
effect of a large number of scalar fields leads to enough number of c-folds. Similar 
idea was exploited in N-flation [6] and Cascade inflation [7, 8] to obtain inflation 
from several potentials that are individually too steep to sustain inflation. In these 
models, although the "effective" inffaton fleld gets a super-Planckian fleld value for 
chaotic m^0^ and inflationary scenarios, due to the large number of flelds, each 
physical fleld remains sub-Plankian. Moreover, for the case of the inflationary 
theory this can be used to resolve problem of unnaturally small coupling [9, 10]. Also 
following the Lyth bound [11], due to the large excursions of the effective inflaton in 
comparison with Mp, in these models one expects to obtain a considerable amount 
of gravity waves. 

In this work we take a different view and promote the inffaton fields to general 
N X N hermitian matrices and hence these models will be called Matrix Infiation or 
M-flation for short. In this sense M-flation is a special case of multi-fleld inflation. 
Working with matrices, besides the simple products of the flelds, we can also consider 
commutators of matrices. In our class of M-flation models we consider three N x N 
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matrices, $j, i = 1,2,3 and the potential is taken to be quadratic in the or in 
their commutator Therefore, in the class of models we consider the potential 

term for $j can have three types of terms: Tr Tr e^jfe^j and Tr <l>^. 

As we will argue, this class of potentials is well-motivated from string theory and 
brane dynamics. 

As we will see despite the simple form of the potential constructed from these 
matrices and their commutators our model has a rich dynamics. We argue that 
M-flation can solve the fine-tunings associated with standard chaotic inflationary 
scenarios. Furthermore, like any multi-field inflation model, there would be iso- 
curvature perturbations as well as the usual adiabatic perturbations. This can have 
significant observational consequences for the CMB observations [1,12]. Moreover, 
we argue that our model has an embedded efficient preheating mechanism. 

The outline of the paper is as follows. In section 2, we provide the setup through 
introducing the action and show that, with the appropriate initial conditions, for the 
sector in which the $i fields fall into irreducible N x N representations of SU{2), the 
theory effectively and at the classical level, behaves like a single field infiation. In 
section 3, we discuss simple models of infiation such as chaotic, symmetry breaking 
and inflection point inflation which are constructed from our M-flation. In section 
4, we work out the mass spectrum of the remaining 3A^^ — 1 scalar flelds, the iso- 
curvature modes, which are not classically turned on during inflation (these flelds will 
be generically called ^f-modes). In sections 5 and 6 we consider quantum mechanical 
excitations of ^-modes and their effects. In section 5, we compute the power spec- 
trum of fluctuations of the infiaton and the 3A^^ — 1 iso-curvature modes. In section 
6, we focus on quantum mechanical creation of the sub-Hubble modes. This particle 
creation takes energy away from the infiaton field. We show that our model naturally 
contains this mechanism [13] which is the basis of the preheating scenarios [14,15]. 
In section 7, we give string theory motivations behind the M-fiation models, arising 
from dynamics of multiple D3-branes in specific fiux compactifications. The last sec- 
tion is devoted to discussion and outlook. In the Appendix we have gathered some 
technical details of slow-roll infiation. 



2. M-flation scenario, the setup 

As explained before, we start with an infiationary setup in which the infiaton fields 
are taken to be iV x non-commutative hermitian matrices. We start from the 
following action 

j d4^^||^i?-l^Tr(a^<I.,a^<I>,)-l-($,, [$,,$,]) j , (2.1) 

where the reduced Planck mass is Mp^ = SttG with G being the Newton constant 
and the signature of the metric is (— , -|-, -|-, -|-). Note that the i-index counts the 
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number of matrices and is not denoting the matrix elements; the matrix element 
indices are suppressed here. Also, V represents our potential constructed from the 
N X N matrices $i and their commutators The kinetic energy for <l>j has 

the standard form and it is assumed that the $i matrices are minimally coupled to 
gravity. 

As we will discuss in section 7, the potential V{^i, can be motivated 

from dynamics of branes in string theory where up to leading terms in $j and , 
in specific string theory backgrounds, the potential takes the form 

V = Tr ^j] + jej,i[^,, + ^^{^ , (2.2) 

where i = 1,2,3 and hence we are dealing with 3A^^ real scalar fields. The Tr is 
over N X N matrices, and here and below the summation over repeated i,j indices is 
assumed. A is a dimensionless number while k and m are constants with dimensions 
of mass. We take A, k, and to be positive. Note that the potential (2.2) is 
quadratic in powers of and The action (2.1) together with the potential 

(2.2) are invariant under U{N) (acting on the matrices) and SU{2) (acting on i,j 
indices) which are both global symmetries. 
Starting with an FRW background 

the equation of motions are 

+ SH^i + A [^i, ]+iK eijki^j, + m'^^i = , 
H = — Y^ Tid $ 9^$ 

^ i 

where H = a/ais the Hubble expansion rate. 

2.1 Truncation to the SU{2) sector 

As discussed, $i are N x N matrices and we are hence deahng with 3A'"^ real scalar 
fields which are generically coupled to each other. This makes the analysis of the 
model in the most general form very difficult, if not impossible. Noting the specific 
form of the the potential (2.2) and that i,j indices are running from 1 to 3, there is 
the possibility of consistently restricting the classical dynamics to a sector in which 
we are effectively dealing with a single scalar field. This sector, which will be called 
the SU (2) sector, is obtained for matrix configurations of the form 

$, = 0(t)J, , i = l,2,3, (2.5) 
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(2.3) 



(2.4a) 
(2.4b) 
(2.4c) 



where Jj are the basis for the dimensional irreducible representation of the SU (2) 
algebra 



N 

[Ji, Jj\ =ieijkJk , Tr(Jj Jj) = ^(^^ " 1) % • 

Since both $j and Jj are hermitian, we conclude that is a real scalar field. 
Plugging these into the action (2.1), we obtain 



(2.6) 



S — d^x\f^^ 



Ml 



-R + Tr 



(2.7) 



where Tr = Y.]=i Tr( Jf ) = N{N^ - l)/4. 

Interestingly enough, this represents the action of chaotic inflationary models 
with a non-standard kinetic energy. Upon the field redefinition 



(Trj2) 



2\-l/2 



-1/2 



(2.8) 



the kinetic energy for the new field becomes standard, while the potential for it 
becomes 



where 



^eff - 



8A 



2A 



Tr J2 N{N^ - 1) ' 



K 



2k 



(2.9) 



(2.10) 



VTrJ2 ^iV(iV2 - 1) 
2.2 Consistency of the truncation to the SU{2) sector 

The SU{2) sector seems to be a special sector of the M-flation action in which the 
theory becomes tractable and very simple. However, we need to make sure that this 
truncation to the SU{2) sector is indeed consistent with the classical dynamics of 
the model and that we can consistently turn off the other SA^^ — 1 fields. In order 
to do this, we define 

*i = - 4>Ji (2.11) 

where, as before. 



= 



:Tr($,JO , 



N(N^ - 1) 

and hence Tr(^j Jj) = 0. In other words, is defined such that it has no components 
along Jj. Using [Jj, Jj] — ieijkJk, we may rewrite the potential (2.2) in terms of (p 
and as 

y = l/o(<^) + l^(2)(0,*i) (2.12) 
where Vq, after the field re-definition (2.8), is given by (2.9). 
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Using (2.6), (2.11) and the fact that Tr(^iJj) — 0, one can show that V(2) does 
not have any hnear terms in i.e. 

y(2)(0,^, = O)=O, =0. (2.13) 

This leads to the important result that the field does not source the fields. 
Explicitly, in the equations of motion (2.4b) if we start with the initial conditions 

= 0, = and (6 ^ 0, will always remain zero and will hence not contribute 
to the classical background inflationary dynamics at all. This means that we are 
consistent in considering as the sole fleld driving the inflation. The remaining 

modes, however, as we shall see in next sections, will be excited at the level of 
perturbations. So, we have SA?"^ — 1 ^ modes besides the background inflationary 
field. In the language of [3] this corresponds to a straight inflationary trajectory in 
the fleld space of 3A'"^ scalar flelds, where the inflationary trajectory is along the 
direction and there would be 3A'"^ — 1 iso-curvature perturbations perpendicular to 
the adiabatic trajectory. We should stress that this result will not remain valid if 
initially fields are also turned on. 

One may wonder about the special role of the S\J (2) generators Jj among other 
N xN matrices in our analysis and whether a similar reduction to a sector other than 
the SU{2) sector is also possible. Traces of the fact that SU{2) sector is special is 
already built in the initial construction of the potential (2.2) with three i = 1, 2, 3 
(3 is the dimension of the SU (2) algebra) and that in the cubic term in the action the 
structure constant of SU{2) Cijk appears. To see this more explicitly, let us consider 
the more generic decomposition for $j matrices 

^^ = T, + E,, (2.14) 

such that Tr(TjSj) = 0, we take both Tj and Sj to be hermitian. The potential will 
again have two parts, V = Vo(Ti) + V(i)(T, 5), where Vo(Tj) = V{Ei = 0) and 



l^(i)(T„S,) = Tr 



-A[Ti, [Tj, Tfe]] + meijfc[Ti, Tj]^ Ek 



+ 0(S2) 



In order the T-sector to decouple, the above expression in the bracket should vanish 
for any Sj. As explained above, if Tj is proportional to Jj this condition obviously 
holds. In general, however, this can happen if the T-dependent term in parenthesis 
is proportional to T^. This condition can be satisfied if [Tj, Tj] = fijk'^k for some 
functions fijk which means that the three Tj matrices should form a Lie-algebra of 
dimension three. The only non-trivial solution is then Tj forming an SU (2) algebra.^ 
Note, however, that the representation for the SU{2) algebra is not fixed by this 
requirement and Tj could also form reducible N x N representation of SU{2). 

^There is also the trivial choice of /^fc = 0, corresponding to choosing three Abelian subgroups of 
U{N). Working with commuting matrices, however, kills all the interesting inflationary dynamics. 
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Motivated by this unique property of the SU (2) sector, from now on we assume 
that the background inflationary trajectory is along the direction, while the other 
3A^^ — 1 fields are only excited at the level of (quantum) perturbations. 

3. Various inflation models resulting from M-flation scenario 

The SU{2) sector, which is governed by the potential (2.9), depending on the values 
of the parameters Ag//, Heff and m^, leads to different inflationary models. These 
models have been studied in the literature, which we will review below. The impor- 
tant advantage of our M-flation scenario, as we will show, is the scaling properties 
between the original parameters of the model appearing in the action and the effective 
dressed parameters appearing in (2.9), which allow us to remove the unnaturalness 
and fine-tuning of these parameters. 

3.1 Chaotic inflation 

If we set A = K = 0, we obtain the simple quadratic chaotic potential and in this case 
our model is nothing but an N-flation model [6] with 3iV^ fields. To fit the CMB 
observations and obtain right number of e-foldings, one needs that m ~ 10^^ GeV 
and a super-Planckian field variation for effective infiaton field during infiation, 
A0 ~ lOMp. In the context of effective field theory this may sound problematic. 
However, as in the N-fiation model [6], note that is our physical field and A0 ~ 
jy--3/2 ^0 -poY a sufficiently large value of one can arrange that A(f) <^ Mp and 
one can avoid the problem with super-Planckian field values. 

On the other hand, if m = k = 0, we obtain the quartic Ae//0'^/4 chaotic 
inflation. To flt the COBE normalization and obtain right number of e-folds, one 
requires Ae// ~ 10"^'' and A0 > lOMp. In the context of a single scalar field, these 
are viewed as severe fine-tunings in the model. In our case, however, we see that 
both of these can be relaxed. Assuming that A ~ 1 dictated by the naturalness 
of the theory, to obtain the above value for Ae// one needs to have ~ 10^. In 
the context of string theory studied in section 7, where N is viewed as the number 
of coincident branes, this is easily achieved in light of recent developments in the 
flux compactiflcation. With this value of N for the physical fleld one obtains 
Act) ~ 10''^ Mp and can hence safely bypass the problem of super-Planckian excursion 
of the fleld. In both of these examples, inflation ends when 0(t) —>■ which means 
that the matrices $j commute with each other. 

In both examples, and for next two examples below, due to super-Planckian 
value of A0 during infiation, a considerable amount of gravity waves can be produced 
which can be detected in future gravity wave observations such as PLANCK [16,17], 
CMBPOL [18] and QUIET [19]. 
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3.2 Symmetry breciking inflation 

Now consider the general case where none of the coefficients in Vo(0) is zero. We 
study two interesting example here. The first example is when Vo(0) is positive 
definite and has two degenerate minima. This corresponds to «; = 3m\/A/2 and the 
potential has the form 

Vo^^ct>'{<t>-Hr (3.1) 

where /x = \^m/ \/Kff- As mentioned, the potential has global minima at = 
and (f) — /I- In the language of the brane construction (c/. section 7), the minimum 
at (f) — II corresponds to super-symmetric vacua when N D3-branes blow up into a 
giant D5-brane in the presence of background RR field C^^^ which in our notation 
is represented by k. The minimum at = 0, on the other hand, corresponds to 
the trivial solution when matrices become commutative. If we allow the field (p 
to take negative values, then the potential is symmetric under (f) —(f) + fi and 
it represents the standard double well potential, justifying the name "symmetry 
breaking" inflation. 

Potential (3.1) is well studied in the literature and is in the form of symmetry 
breaking potential and "hilltop inflation" [20,21]. In Appendix A we briefly look 
at the predictions from this potential compared to WMAP5 data. It is assumed 
that inflation lasts for 60 number of e-folds, A'^e = 60, and the scalar spectral index. 
Us = 0.96, from WMAP5 central value. Furthermore, the COBE normalization is 
set to 5h ^ 2.41 x 10"^ 

Depending on the initial value of the inflaton field, 0j, the infiationary period is 
divided in three categories. 

(a) (f)i> n 

Suppose inflation starts when 0j > fj,. With Ng — 60, Sh — 2.41 x 10~^ and Ug — 0.96, 
one obtains 

0i~ 43.57Mp , 0/~27.O7Mp , ~ 26Mp. (3.2) 

and 

Xeff ~ 4.91 X 10-^^ m ~ 4.07 x 10-^Mp, Keff ~ 9.57 x 10-^^Mp. (3.3) 

(b) /i/2 <(t)i<ii 

This is an example of hill-top inflation when the infiaton field is between the local 
maximum at = /x/2 and the "supersymmetric minimum" at = /x. To fit the 
above observational constraints one obtains 

(pi ~ 23.5Mp , (t)f ~ 35.03Mp , ~ 36Mp. (3.4) 

and 

Ae// ~ 7.18 X 10-^^ m ~ 6.82 x lO^^Mp, K^ff ~ 1.94 x lO'^^Mp. (3.5) 
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(c) Q<(i)i<^l/2 

Due to symmetry (f) —0 + fi this inflationary region has the same properties as 
11/2 < (bi < ji above. If we allow for negative values of 0, then the inflationary 
prediction for </> < region is the same as in region 3.2 above. 

In all these examples, to fit the COBE normalization, one obtains N ~ 10^ 
and ~ 10~^Mp during inflation so the issues with super-Planckian field range is 
resolved. 



3.3 Inflection point inflation 

Another class of models widely studied in the hterature coming from the potential 
(2.9) is when the potential has an inflection point [22]. This happens when k — 



v2Am or equivalently, Kg// — m-y/Ag//. Denoting the inflection point by the 
potential near 0o is approximately given by 

V{(t>) ~ V^(0o) + ^V"\(t>o){(t> - cPof (3.6) 



where 



n0o) = ^0^ , ^'"(0o) = ^. (3.7) 
12 00 

The CMB observables are given by [22] 

= , fe^AMMfljv^^ (3.8) 

Choosing = 60, one obtains ~ 0.93 which is within 2a error bar of WMAP5 
but is somewhat to its lower end. One may add small modification to the coefficient of 
potential [22] such that the potential around the inflection point is slightly modified. 
This in turn can result in a higher value of Ug. On the other hand, from COBE 
normalization, one obtains 

V/ ~ 10"' ■ (3.9) 
" Mp ^ ' 

In a conservative limit that m < Mp, this yields Ag// < 10~^. Starting with A ~ 1, 
this corresponds to A^ > 10^. 

If we look into the amplitude of the gravity wave, determined by quantity r which 
is the ratio of gravitational perturbation amplitude to scalar perturbation amplitude, 
one obtains 

. fV'V 2 ( m 
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Combining this with Ug and 5hi we obtain 

Ae// = (^)^/^(f<^^)'(l-n.r/^ (3.11) 

The upper bound, r < 0.22, from WMAP5 imphes that Ae// < 10"^^ and N > 10^ 
which is stronger than the bound above. 



4. Mass spectrum of modes in M-flation 

Having studied the SU (2) sector and the resulting inflationary models, we review our 
model, noting that the other 3N'^ — 1 fields encoded in ^j, although not contributing 
to the classical inflationary dynamics, do have quantum fluctuations and will hence 
affect the cosmological perturbation analysis. To compute these effects we need to 
have the mass spectrum of the 3A^^ — 1 modes coming from the ^j. 

To this end, starting from (2.11) we expand the action up to the second order 
in ^'j. Given the orthogonality condition, Tr(\E'jJi) = 0, the kinetic term readily 
takes the standard form |Tr(9^^i9^\E'i). After a shghtly lengthy but straightforward 
computation the potential to second order in is obtained as 



V^2) = Tr 

where 



(4.1) 



nk = ieijk[Ji,'^j] ■ (4.2) 



^Prom the above form we see that if we have the eigenvectors (eigen-matrices) 
of the Qj we can compute the spectrum of in terms of 0-field (to be viewed as 
the inflaton). Finding the eigenvectors of fli is mathematically the same problem 
as finding the vector spherical harmonics. (For a detailed discussion see e.g. [23], 
section 5.2.) If we denote the fl eigenvalues by lu, i.e. 

= ou^i , (4.3) 

we obtain 

If we have the possible values of uj we can read off the effective (0-dependent) mass 
of the modes 



V^'iuj + if -^{4uj + 3){uj + 2) + -J-{uj + l)(a; + 2) , 
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where Vq and V^q" denote the first and second derivatives of the potential Vq with 
respect to the infiaton (f). For later convenience it is also useful to write the expression 
for the mass during inflation in terms of the Hubble expansion rate and the slow-roll 
parameters e and rj 

rj{u; + 1)2 - sgn{V^)V2'€ ^{Au + 3)(a; + 2) + 6^(a; + l)(a; + 2) 

0^ 



(4.6) 

where sgniyl) represents the sign of Vq and as usual the slow-roll parameters are 
defined by 

Following the analysis of [23] , we find that to can take three values: 
• " The zero modes" uj = —1. This happens for modes of the form 

with A being an arbitrary matrix. For these modes the expression for the mass 
simplifies to 

= ^ . (4.8) 

Therefore, at the minimum values for where Vq vanishes these modes become 
massless. (This justifies the name "massless modes".) 

Noting that A is an arbitrary matrix there are N"^ of such modes, all with the 
same mass. 



"The a modes": u ^ -(I + 1), I E Z, 0<l < N, with the mass 

2 



Mf = ^{l + + 2)02 - 2Keff{l + 1)0 + m2 . (4.9) 



Each mode for a given I has a multiplicity of 2Z -|- 1. Therefore, there are N'^ 
of CK-modes. 

• "The f3 modes": u ^l, I e1^, < I < N, with the mass 

= ^/(Z - l)02 + 2«;e///0 + m2. (4.10) 

Each mode for a given I has a multiplicity of 2Z-|- 1. Therefore, there are — 1 
of /3-modes. 

As expected, there are altogether SN"^ — 1 zero, a, and (3 modes. For the zero 
modes where u; — —1, one observes that M'^/SH'^ ~ ^/eMp/cf) ~ 0.01. So there are 

of such light zero modes. For a and /3 modes, only those with I €' 
are light, while the modes with higher values of I are heavy. 
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5. Power spectra in the presence of modes 

With the mass spectrum for modes computed in the previous section we can com- 
pute the power spectra of the adiabatic and the iso-curvature perturbations. Here 
wc arc deahng with a 3iV^ real scalar field inflationary system. Our inflationary 
background is along the direction. The remaining 3iV^ — 1 scalars are frozen classi- 
cally during inflation and are excited only quantum mechanically. Correspondingly, 
the modes are classified as the adiabatic perturbation, the one which is tangential 
to the classical inflationary background, and the iso-curvature modes, the 3N^ — 1 
perturbations which are perpendicular to the background inflationary trajectory. In 
this respect our model is similar to the model studied in [24] where the potential 
has an 0{N) symmetry such that the inflaton field is the radial direction while the 
remaining N — 1 angular directions are iso-curvature perturbations. 

The formalism to calculate the adiabatic and iso-curvature (entropy) power spec- 
tra was systematically developed in [3]. Here we shall repeat those analysis for our 
system of 3A'"^ scalar fields. 

5.1 Lineeir perturbations 

The perturbed metric in the longitudinal gauge is 



where $(t, x) is the gravitational potential and should not be mistaken with which 
represents our matrix fields. Similarly, the linear perturbed scalar fields are 6(j){t, x) 
and ^i(t,x). 

We find it is much easier to work with the modes with definite mass spectrum. 
As discussed in the previous section, these fall into three classes denoted by ^r,im 
where r — 0,a, P stand, respectively, for zero-mode, ct-mode and /3-mode and m is 

running from one to Dj.j, the degeneracy factor of each '^r,im. mode, with Dq = iV^ 
and Da,i = = 2/ + 1. In this notation, the Lagrangian for the scalar fields with 
potential V = Vb(0) + V(2)(0, *i) from (4.4) is 



ds'^ = -(1 + 2<^)dt^ + a{tY{l - 2<^)dx^ 



(5.1) 



(5.2) 




(5.3) 




(5.4) 
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Here k is the momentum number in the Fourier space. One interesting aspect of 
these equations is that the adiabatic and iso-curvature modes completely decouple 
and they do not source each other. This is a result of our initial conditions in which 
"ifrM fields are turned off. Because of this, the equation of motion for ^r,Zm is that 
of a scalar field in a homogeneous and isotropic expanding background but with a 
time dependent mass. 

These equations are accompanied by Einstein equations 

^ + ^$=^050- (5.5) 

Interestingly enough, ^r,im, fields do not show up in perturbed Einstein equations 
due to our assumption that they are absent at the background dynamics. Physically, 
this means that the modes ^r,im do not carry energy up to linear order in pertur- 
bation theory. Since they do not couple to gravitational potential and the inflaton 
perturbations, they have no gravitational effect, justifying the name iso-curvature 
perturbations. 

Consider the normalized curvature perturbation TZ and the normalized iso-curvature 
perturbations 5r,/rn 

H H 

TZ = -j-Qfli , Sr,lm = -;-^r,im • (5.6) 
(f) (f) 



Using the Einstein equations one can show 

7t=^V- (5-7) 

This indicates that on arbitrary large scales where /c — > 0, the curvature perturbation 
is conserved. This is similar to single field inflationary system. 

Here we assume that initially Scf) and ^r,;m ^-re random, Gaussian and adiabatic 
fields which are excited quantum mechanically from vacuum and 

(Q;kg<^k') = ^^Q/(k-k') 

27r2 3 / 

{^r,lm k *r',/'m' k') — -j^P^^^i ^rr' ^ll'^mm' ^ (k — k ) 

(g^Aimk') = 0, (5.8) 

where Pq^ and P^^ ^ are the primordial power spectra of the scalar fields. We note 
that the last equation above also holds during the inflationary stages, indicating that 
there is no cross-correlation between adiabatic and iso-curvature modes. Physically, 
this means that they do not source each other as can be seen from (5.4). 
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5.2 Power spectra from Hubble-crossing to end of inflation 

As usual, it is instructive to express the Klein- Gordon equations (5.4) in terms of 
variables u and Vrm defined by 



u = aL 



Vr,lm = a^r,, 



Im 



Going to conformal time dt — adr, (5.4) is transformed into Schrodinger equation 
form 



dr^ 



+ 



2 - 377 + 9e 



u = Q 



d?u 



r.lm 



dr^ 



+ 



2 2 - 377^,/ + 3e 



VrJm = 



(5.9a) 
(5.9b) 



where r]r,i = M!^i{(f))/3H'^ given by (4.6). 

The initial conditions deep inside the Hubble radius are given by the Bunch- 
Davies vacua. We evolve equations in (5.9) till the time of horizon crossing r* at 
which k — (aH)^. Here and below the subscript * indicates that the quantities are 
calculated at the time of Hubble crossing during inflation. Following the standard 
inflationary prescriptions e.g. [25], up to the first order in slow-roll parameters one 
has 



u 



TT r 



ik\r\) 



(5.10) 



where H^^^ (x) is the Hankel function of the first kind and ur and vs^ are given by 



(5.11) 



Correspondingly, using (5.6), the power spectra of the scalar curvature perturbations 
and iso-curvature perturbations, P-ji and Ps^ ,^ , at the time of Hubble exit are given 
by 



111* 



H 



2 \ 2 



27r0 
2 \ 2 



[1 + (-2 + 6C)e - 2Cr]l 



Psr.J* ^ [1 + (-2 + 2C)e - 2C7ir,i] 



(5.12) 



Here C = r'(3/2)/r(3/2) + ln2 ~ 0.7296 where r(x) is the Gamma function. 

Few e-folds after the mode of interest has left the Hubble radius, one can neglect 
the term k'^/a^ in (5.4). Furthermore, it would be more instructive to use the number 
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of e-folds before the end of inflation, N^, as the clock. Using the relation dN^ — Hdt, 
(5.4) is transformed into 

^c^{2e-r))dN, , ^c^-r^,,dN,. (5.13) 

To calculate the power spectra as a function of Ng, for TZ and Sr^im defined as in 
(5.6), we need to incorporate the evolution of the pre- factor H/(j) which is given by 

' ^''^^-{2e-rj)(^). (5.14) 



dN, 

Combining equations (5.13) and (5.14) one obtains 

Pn{Ne) ~ PnU (5-15) 



-2 / dN'^Br,i{K] 
Jo 



(5.16) 



where Br,i{Ne) =2e — rj + rjr^i and using (4.6) one obtains 



B,,i{N,) ~ 2e + o;(2 + a;)77 - sgn{V^)^^^^{Auj + 3)(a; + 2) + ^(a; + l)(a; + 2) 

where uj takes values 0, —l — l and I respectively for zero, a and (5 modes (c/. section 
4). 

As explained previously, due to conservation of IZ on super-Hubble scales, P-ji 
remains unchanged after Hubble exit as shown in (5.15). On the other hand, the 
evolution of power spectra for the iso-curvature modes, Ps^i^{Nf.)^ as shown in (5.16), 
depends on the dynamics of the inflationary background and the eigenvalues uj. 

5.3 Iso-curvature vs. entropy perturbations 

As we saw in the previous section, the '^r,im perturbations do not couple to the 
inflaton field and the gravitational potential so they do not contribute to the primor- 
dial curvature perturbations during inflation which justify the name iso-curvature 
for these modes. Physically, this means that up to the first order in perturbation 
theory, ^r.im fields do not carry energy during infiation. As showed in section 2.2 
this has the origin in our initial conditions for the '^r,im fields that they are absent 
in classical background dynamics. 

In the literature, the terminologies "iso-curvature perturbations" and "entropy 
perturbations" are usually used interchangeably. However, one can easily check that 
^r,im perturbations do not induce entropy perturbations during inflation. To see this 
explicitly, let us look at non-adiabatic components of pressure, Spnad, defined as [4] 

Spnad^ Sp-^Sp, (5.17) 
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where p and p are the background energy density and pressure respectively and 6p 
and 5p represent their first order variations. Using Einstein equations (5.5) one can 
show that 

5p„„,~-4M|^$, (5.18) 

which is the same as the standard single-field infiation result. On super-Hubble scale 
one observes that Spnad vanishes and there is no entropy perturbation. This is similar 
to our earlier result (5.7) that 7?. ~ on super-Hubble scale. In this work, in order to 
emphasis that there is no entropy perturbations in our setup, we use the terminology 
iso-curvature perturbations throughout. 

In general multiple-field infiation, the iso-curvature perturbations perpendicular 
to the classical infiationary trajectory do produce entropy perturbations and Spnad 7^ 
0. These entropy perturbations source the adiabatic perturbation and have non-zro 
energy at the linear order of perturbations so they also source the Einstein equations 
and contribute to the primordial curvature perturbations. More specifically 

7^ = ^ — $ + 2 V 9r,lmSrM ■ (5.19) 
r, Im 

The above equation is a generahzation of the two-field result of [3] , the sum in the last 
term is over — 1 entropy modes and 9r,im represent the angle between ^r,/m and 
the infiaton trajectory, cf) in our case, in the space of 3A^^ scalar fields. As discussed 
above, if we start with the initial conditions \E'j = 0, = 0, behaves as the sole 
infiaton field and in the language of [3], the background infiationary trajectory in 
— phase space is fiat, i.e. 0^,;™ = 0. Therefore, these iso-curvature perturbations 
do not feed the curvature perturbation and there is no cross-correlation between the 
iso-curvature and adiabatic perturbations. 

However, if we start with an arbitrary initial condition in ^ — \l'j field space, ^'j 
are not frozen classically, the infiation trajectory is curved and the infiaton field has 
a component along the "ifrM direction and 9r,im 7^ 0. This in turn produces entropy 
perturbations during infiation which also source the cosmic perturbations. 

Having this said, however, there are two mechanisms to create entropy pertur- 
bations in our setup. The first way is to consider second order perturbation in 
^r,Zm- In second order perturbation theory, ^r,Zm carries energy [27] and couple to 
both Einstein equations and the infiaton field perturbation 6(f). This in turn leads 
to entropy perturbations during infiation. However, the amplitude of these entropy 
perturbations are much smaller than the first order adiabatic perturbations coming 
from 50. The second mechanism to create entropy perturbations from iso-curvature 
perturbations "ifrM can happen during preheating and/or reheating era. Similar idea 
was used in [24] through asymmetric preheating. 
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Figure 1: Left and right graphs respectively show the curvature and isocurvature spectra 
for chaotic inflation with potential 

5.4 Curvature and iso-curvature perturbations power spectra for specific 
examples 

Using the general formulation presented in previous subsections we calculate the 
power spectra at the end of inflation for chaotic, symmetry breaking and the inflection 
point inflationary potentials. 

5.4.1 Chaotic infiation 

"When A = AC = 0, the action takes the form 

S = ~d,<pd^<f> - ld,%^^d^%,ir, - [0^ + K,J^rM] ■ (5-20) 

The 3A^^ modes all have equal masses. Interestingly, this closely resembles the chaotic 
infiation studied in [24,26]. The potential (5.20) has 0(3A^^) symmetry for the fields 
(f) and ^r,;m- The background inflation field is and the number of e-folds in terms 
of the initial value of the inflaton field 0j is given by 



4M|. ■ 



Correspondingly, the amplitude of the adiabatic curvature perturbation is 

2 



1 



mNr- 



67r2 V Mp 



(5.21) 



(5.22) 



For Ne — 60, to fit the WMAP5 normalization P-ji ~ 2.41 x 10~^, one requires 
m ~ 6.304 X 10~^Mp. For such values of m, we have calculated numerically the 
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Figure 2: Left graph shows the ratio Ps^i^i^e)/ Ptz\* vs. N^, for the mode that exit the 
Hubble radius 60 e-folds before the end of chaotic inflation with V = ^m?(lP'. Right graph 
shows the evolution of Psf^ ^^{N(,) / P-jiU for Z = 1 /3— mode in Ae//0^/4 potential. Black 
and gray curves respectively demonstrate the numerical and analytic results. 

amplitudes of 3A^^ — 1 iso-curvature spectra at the end of inflation, i.e. where e = 1, 
which turns out to be Psrim — 8-426 x 10~^^ at today's Hubble scale. The spec- 
tral indices of adiabatic and iso-curvature spectra at such scale, respectively, are 
tiTi ~ 0.966 and n^^^^ ~ 0.9998, which coincide with the analytic results of [24], see 
Fig.l. One can lower the value of m - which in turn lowers the amplitude of adia- 
batic perturbations - and generate the difference by transforming the iso-curvature 
fluctuations to curvature ones through an asymmetric mechanism of preheating [24]. 
The ratio Ps , I P-n for the mode that exit the horizon 60 e-folds before the end of 
inflation is graphed as a function of A^e, see Fig. 2. Note that as in this case all the 
^r,im modes have the same mass and hence 77^,;, Psri^jPn independent of r and 
Z, m. When the mode is inside the Hubble radius the spectra are almost equal and 
the ratio is one. However, around the Hubble-exit, the iso-curvature mode starts to 
decay following (5.16). For the mode that exits the Hubble radius 60 e-folds before 
the end of inflation, /cgo = e~^°aei?e, according to (5.16) the ratio decays like 

^g" ' ~ (1 - A^e/60)' , (5.23) 

where N^. is the number of e-folds the mode spends outside the Hubble radius before 
the end of inflation. As can be seen in Fig. 2, the analytic result is in a good 
agreement with the numerical one. 

For this model, the amphtude of tensor fluctuations of the adiabatic mode at 
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today's Hubble scale is Pri^m) — 3.1856 x 10^^°. The corresponding tensor/scalar 
ratio, r, isr ~ 0.132. The tensor spectral index, rir = din Pt/ dink, isnx — —0.0165. 
Future CMB probes such as PLANCK [16], or exclusive polarization probes such as 
CMBPOL [18] or QUIET [19] should be able to test this scenario. 

5.4.2 Chaotic inflation 

When m — K, — 0, the potential energy V — Vo{(f)) + V(2)(0, ^i) is 

V = ^<P' + ^(o;^ - c.) KM'^rM ■ (5-24) 

To match the amplitude of adiabatic perturbations with WMAP result at horizon 
scale, Xeff — 1.6315 x 10~^^. For such value of Xeff, the scalar spectral index for the 
adiabatic mode is n-R, ~ 0.949. 

In Xeff(p'^/4: chaotic model, the masses of iso-curvature modes are different. The 
lowest mass belongs to the I — 1 /3— mode whose mass is equal to zero, M|^(0) = 0. 
The corresponding iso-curvature spectrum amplitude at the end of inflation is equal 
to Psg 1^ — 3.949 X 10~^^ at today's Hubble scale. The corresponding spectral index 
is riqig ~ 0.966. The relatively larger value of iso-curvature perturbations could be 
attributed to the fact that the iso-curvature spectrum for this mode decays linearly 
with the number of e-folds it spends outside the horizon (sec below). 

The next modes in the tower of masses are the zero mode, / = a— mode and 
I — 2 /3— mode whose masses are equal to Mg 2(0) = ^eSf4>^- Their corresponding 
amphtudes are Psg ^ 4.449 x 10~^^ at today's horizon scale. Taking the multi- 
plicity of a and (3 modes into account, there are N'^ + 6 iso-curvature modes with 
such an amphtude. The corresponding spectral index is ~ 0.9828. 

For the / = 1 a— mode with the mass M|3(0) = 3Ae//0^, which is equal to 
the mass of the / = 3 /5— mode, the iso-curvature spectrum amplitude is Pg^ .^^ — 
3.967 X 10^^^ at today's Hubble scale. The corresponding spectral index is equal to 
n^^jj^ ~ 1.016, which indicates a blue spectrum. 

In general mass of a / > 1 a— mode is identical to the / -|- 2 /5— mode. Therefore, 
there are 4Z -|- 6 iso-curvature modes with identical spectra, all of which have a blue 
tilt. Increasing the value of I for ct— mode and /3— mode, the amplitude of heavy 
iso-curvature modes decreases quickly. This can be understood from (5.16), which 
yield 

iVe/60)^ zero modes 

^^/gQ^(z^+3Z+4)/2 «_niodes (5.25) 

Nj60f'-^^^^^^ (5 - modes, 

where again is the number of e-folds the mode spends outside the Hubble radius 
before inflation ends. In particular, for Z = 1 /3-mode our analytic result indicates 



^■^\ ^ ~ (1 - iVe/60)'+ — 

Pn\* 
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that the ratio Ps^ ^^{N^) / P-r\* decreases hnearly with N^, which is verified numeri- 
cally as shown in Fig. 2. 

Besides a considerable iso-curvature/adiabatic ratio of 1.638% for / = 1 /9— mode, 
another signature of the model is its observable gravity waves. For this model the 
amplitude of tensor spectrum for adiabatic perturbations at Hubble scale today is 
Pt(A;o) ^ 6.3176 x 10"^°, i.e. r ~ 0.26, whose spectral index is ut ^ -0.033. This 
model is currently on the verge of becoming ruled out. 

5.4.3 Symmetry breaking inflation 

We consider (f) > /jl and ///2 < < cases separately. 

(a) (j)> II 

In this case neither of the parameters of the potential. A, k and m are zero. To 
match the observational constraints from WMAP5, the above parameters have to 
take the following values given in (3.3): 

Ae// ~ 4.91 X 10"l^ m ~ 4.074 x lO'^^Mp, n^ff ~ 9.574 x lO^^^Mp. 

The scalar spectral index for the adiabatic perturbations is n-jz ~ 0.959. 

The lowest masses in the tower of ^r,Zm iso-curvature modes belong to the zero 
mode and the / = «— mode, whose mass is M^o(0) — '^eff4>^ ~ '^'^eff4' + The 
amplitudes and spectral indices of these two iso-curvature spectra are respectively 
Ps „ ~ 1.162 X IQ-ii and „ ~ 0.981. 

Next in the tower of iso-curvature modes is the / = 1 a— mode whose mass is 
equal to ^{(j)) = 3Ae//0^ — 4fi;0 -|- rri^. Its amplitude and index are, respectively, 
Ps^arr. - 6-966 x lO-^^ and n^^,,^ ~ 1.01. 

The I — 1 /3— mode with mass equal to 2Keff(t> + stands next in the tower. 
The amplitude of this mode is equal to Psg ^ 8.842 x 10~^^. The corresponding 
iso-curvature spectrum for this mode has a blue tilt but an almost scale-invariant 
spectrum, with ~ 1.002. As before, the next iso-curvature modes have negli- 

gible amplitudes at Hubble scale and therefore their contributions could be ignored. 

The amplitude of tensor spectrum at Hubble scale is Pt(^6o) — 4.84 x 10~^°, 
r ~ 0.2 with the spectral index ut — —0.025. Planck [16] should be able to verify 
this model. 

(b) ii/2 <(f)< II 

Here to satisfy the constraints from the amplitude and spectral index from 
WMAP5, one has to adjust the the parameters as in (3.5): 

Ae// ~ 7.187 X 10-^^ m ~ 6.824 x 10~^Mp, Keff ~ 1.940 x 10-^^Mp. 
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The index of the adiabatic spectrum for such values of parameters is nji ~ 0.961. 

Again the least massive iso-curvature modes are the zero mode and the / = 
a— mode. Their amplitude and spectral index are, respectively, Ps^ — 1.46 x lO"-*^^ 
and n^^ o^ ~ 0.987. The next biggest iso-curvature amplitude belongs to the / — 1 
a— mode whose amplitude and spectral index are Psc,,im — 9-99 ^ 10~^^ and — 
0.988, respectively. The I — 1 /3— mode stands in the next rank with an amphtude 
of Ps„ , ^ 6.558 X 10~^^ with , ~ 1.054. Next iso-curvature modes have larger 
masses and therefore, as before, their amplitudes are negligible in comparison with 
the adiabatic one. 

The amplitude of tensor spectrum at Hubble scale is Pt(^o) — 1-307 x 10~^^, i.e. 
r ~ 0.048 with the spectral index tit — —0.006. Such gravity wave spectrum could 
be detected by CMBPOL [18] or QUIET [19]. The tensor spectrum is very close to 
being scale-invariant in this case. 

5.4.4 Inflection point inflation 

A possible parameter set that satisfies the observational constraints are: 

Xeff ~ 4.8 X 10-^^, m ~ 10~^Mp, K^ff ^ 1-94 x lO^^^Mp. (5.26) 

For this parameter set the scalar spectral index is n-jz ~ 0.93 which is within 2a error 
bar of WMAP5 but is somewhat to its lower end. 

The lowest mass mode belongs to the / = 1 a— mode whose mass is M^^{(j)) — 
3Xeff(f)^ — 2K.eff(f) + 'm'^. Its corresponding spectrum amplitude and spectral index are 
Ps^^irn. — 1-38 X 10~^° and — 0.932, at Hubble scales which is interestingly very 

close to the spectral index of the adiabatic spectrum. This mode has the highest ratio 
of Psoc,imi-^e) / Pn\* — 5.7%. The next modes in the series are the zero mode and the 
I — a— mode whose corresponding amplitudes and spectral indices are respectively 
Pc ~ 1.7551 X 10^^^ and ^ — 1.0012. Other iso-curvature modes have 
completely negligible spectra in comparison with the adiabatic one. 

The amplitude of the tensor spectral index is negligible for the above set of 
parameters in this model, Pt ^ 1.168 x 10~^^, which is scale-invariant with the 
precision of 10~^. 

6. End of inflation and preheating 

As discussed while the field is turned on during inflation, the other fields ^r,Zm are 
also present. Although not turned on at the onset of inflation (by the choice of initial 
conditions) and hence due to the speciflcs of the classical dynamics of our model 
remain zero during inflation, the presence of '^r,im flelds can be felt through quantum 
effects. These quantum effects show up in two different contexts; one is of course 
through the power spectrum of the quantum fluctuations of these fields at the super- 
horizon scales, the iso-curvature modes which were discussed in some detail in the 
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previous section. The other quantum effect is the possibihty of creation of the \E'r,Zm 
particles, due to the coupling to the inflaton field (j). If the pair creation mechanism 
is "efficient enough" this will eventually back react on the classical dynamics of the 
system. This effect, if too efficient during inffation and before completion of the 
needed 60 e-folds, can tamper the whole M-fiation scenario.^ RecaUing the large 
number of ^r,im modes {?>N'^ — 1), their collective effect on the inflaton field could 
be very large ending inflation too fast. While if activated only toward the end of 
inflation it will be a positive feature of our model, providing us with a mechanism 
to end inflation while transferring the potential energy of the inflaton field into the 
kinetic energy of the ^r,/m fields, a preheating scenario [14, 15,25,29]. 

In the first subsection, we first show that particle creation during slow-roll in- 
flation is not harmful to our M-flation model. In the next subsection we explore the 
possibility of the particle creation as the basis for a natural and inherent preheating 
scenario in our M-flation model. 

6.1 Pctrticle creation during slow-roll inflation is not hctrmful to M-flation 

In this section we show that quantum production of '^r,im modes during inflation is 
not large enough to derail the slow-roll M-fiation. This is done at two steps, first we 
compute the particle creation rate during infiation and then study the back reaction 
of the ^r,;m modes on the dynamics of the inflaton 0. 

6.1.1 Quantum production of '^r,im modes during infiation 

The Lagrangian governing the dynamics of our model, up to the second order in '^r,im 
flelds, and the corresponding equations of motion are given respectively by (5.2) and 
(5.4). As explained before, the time-dependence in Mr^i{4>) leads to '^r,im particle 
creation from the vacuum [13-15,25,29]. (Note that the expression (4.5) for M^i{(f)) 
depends on r and I as well as the momentum number k.) 
By replacing 

Xr,lm k — O^^^ ^r-,im k , (6.1) 

the equation of ^r,im field takes the form of 

Xk + ^k,ri{tf Xk = , (6.2) 
which is an oscillator with a time dependent frequency, fil ri 

Kri = ^'lH^il-le)+M^0), (6.3) 

^This "slow-down" effect via partiele ereation, although potentially harmful for the standard 
slow-roll inflation, can be used as the mechanism to render an otherwise fast-roll inflationary scenario 
which does not give enough e-folds, to an effectively slow-roll inflation with enough number of e- 
folds. The recent publication [28] discusses this possibility. 
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where e = —H/H"^ which in the slow-roll limit reduces to its conventional form (4.7). 

During inflation and when fi^ < and x is either the inflaton or iso-curvature 
perturbations, x modes have an imaginary frequency. As is well-known < 
during inflation (e.g. see discussions and analysis of section 4) happens for the 
"super-horizon" modes. These imaginary frequency modes are those which follow a 
classical dynamics and contribute to the power spectrum of curvature or iso-curvature 
(entropy) modes. 

For particle creation inside the horizon which is what we will be mainly concerned 
with here and can happen if the time variation of VP' is not negligible, we should 
focus on the Q? > regime. Let us suppose that we have the solution to (6.2) for 
the Q'^ > regime. Our goal is to compute the number density of the Xk particles 
produced during inflation and for that we need to compare the solutions at t = — cxo 
(the onset of inflation) to t — +oo (the end of inflation) . Denoting the former by Xk 
and the latter by x^, one may expand 

Xk = ^kXt + BkXl^ ■ (6.4) 

The Ak and coefficients maybe thought as the Bogoliubov transformation pa- 
rameters. The number density of Xk mode produced is then equal to \Bf;\^. For a 
detailed discussion on this matter see [14] or Appendix B of [38] . 

Except for some specific cases (e.g. see [25,30]), however, it is not possible to 
solve (6.2) and we are hence forced to use approximations. One of the approximations 
which is usually employed (e.g. see [15,25,29]) is the stationary phase approximation, 
which if (6.2) viewed as the Schrodinger equation, it is the WKB approximation. This 
approximation is valid if 

^k,rl<^lrl , ^lrl>0. (6.5) 

In this regime one may expand fl^ around its minima as 

Krl = ^krl + ^Irli^ " tkf + 0{{t - hf), Qfe,w|t=t. = . (6.6) 

where tk is where flk,ri vanishes. Note that within our assumptions f^ofc.rz ^k,ri 
are both positive. (If fil ^i several minima one should sum over all of them.) 

In the WKB approximation the number density of the particles produced for 
each mode k is 

iXklXk) = W = e""^ . (6.7) 

As we sec the particle creation is effective if is of order one and not large. 

In order to be specific and to get a better theoretical understanding of the analysis 
we shall focus on the specific Ae//0^/4 model for which Ke// and both vanish. In 
this case 

= uXeffcf)^ , (6.8) 
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where ly — 1, + 1), — 1) respectively for zero, a and (3 modes. Moreover, 



2- 

2 fMp^ ^ 
e — -rj 





In the leading order in e, rj 

1 „o 3 



am^ 21/4^1^-1^(3 ue), ^^^^ 

The creation of particles can then happen in the window 

^ < i/e < 3 . (6.10) 

That is, it is not possible for zero modes and for a. and /? modes with / less or of 
order \l\f^- Recalling the expression (6.7) and the exponential suppression by the 
factor of r^o/^' particle creation is effective in the region 



2 7/^: - 5 



The above can be satisfied if v is close to its lower bound ^. However, recalling that 
V is integer- valued and that for large z/, z/ ^ /^/2, the number of allowed Ts is very 
limited and as a result ve can be tuned around the center value 3/2 with accuracy 
of order ^/e. That is, for v around 

i/e=^ + 5v/^, (6.12) 

with 5 being an order one number, the particle creation is effective. 

As the last step we compute the energy which is carried by the ^r,im particles 
in a and /3 modes produced during inflation. For that we need to integrate over the 
number density (6.7), explicitly 

^ J {27ry a-^Qok,ri 

where sum over r only runs over a and (3 modes and D^^i is the degeneracy of the 
modes which is equal to 2Z-I-1 for a and P modes. Recalling that the particle creation 
is effective in the range, 

a V 8 8 V 2 

and for ue — ^, D^^i is roughly Given the above one can perform the integral 

~ _L ^5/4 ^1/2 . _ /g_;L4) 

We would like to remark that as is of order unity one may still trust the WKB 
approximation. 
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6.1.2 Back reaction on the inflationctry dynamics 

The created \& particles during slow-roll inflation will back react on the dynamics of 
the inflaton 0. Their back reaction can be traced through their effect in the equation 
of motion of the inflaton. Strictly speaking the back reaction effects we want to 
consider arise from the one loop correction to the inflaton potential. Noting the 
action (5.2) and the form of the potential Vq, these corrections are 

4> + 3H^+Vo + AV ^0 , (6.15) 

where ^ 

AV = 3KfM^') - 2Keff{<f>') + - ^.(Klm'^rM)^^ (6-16) 

r,lm 

where (^*;„^r,Zm) during inflation should be replaced with (6.14). The {(f)'^) 
during slow-roll inflation where e, r] are very slowly varying (and are basically con- 
stant), is negligible and can be dropped."^ This term, however, may become large 
toward the end of inflation when both e and e/H can become order one. This is the 
regime we consider in the next subsection. During the slow-roll inflation therefore, 
AV ^iN^dM^ {(/))/ d(l). 

Let us focus on the \eff(t>'^/^ theory for which the computations of 6.1.1 was 
mainly carried out. In this case the back-reaction equation of motion for the inflaton 
0, is 

+ + Ae//0^ + uXeffN^ = 0. (6.17) 

As the computation of A^,^ was carried out assuming the slow-roll approximation 
is valid, we stress that the above equation is hence only trustable during slow-roll 
inflation. 

To check whether the last term is harmful to the (slow-roll) inflationary dynamics 
we compare the last term to contribution of the main potential driving inflation: 

Xeff Z/0 ■ _ 27 f /4 ^1/2^-^5 (a i 



^Note that at the one loop level besides AV terms one should consider the running of the 
coupling constants Xeff, Heff and m?. Since we are dealing with quadratic potentials which are 
renormalizable, one may compute these loop corrections. Due to the large number of fields these 
one loop effects could be large. In our case, however, we can just use the renormalizcd values for 
these parameters. We also note that among the specific cases that we discussed, the Ac/)'' theory at 
one loop level will receive a contribution of the form 5'm?(jP' . The condition of having an inflection 
point, Kgjj = w?\eff will not be preserved by the quantum corrections (at one loop level). The 
symmetry breaking potential, however, will preserve its form. Besides the renormalization of the 
couplings there is also the Coleman- Weinberg corrections. In our discussion, however, we will not 
consider this. The symmetry breaking case can be viewed as the bosonic part of the potential in a 
supersymmetric theory, for which the Coleman- Weinberg correction is absent. 

^The {(fP') is essentially the same {^'^) as the zero ^E* modes for which the v factor, which causes 
the enhancement in the a and (3 modes of is absent. Therefore, the {(jP) is negligible compared 
to {K,lm) term. 
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Recalling the WMAP bound H/Mp < 10"^ the above expression is of course much 
smaller than one during slow-roll inflation. We therefore conclude that the particle 
creation is not going to destroy our M-flation model. 

6.2 Pcirticle creation and the preheat scencirio 

We argued that during the slow-roll inflation particle creation, and hence its back- 
reaction on the dynamics of the field is not large. However, particle creation can 
become important when e, rj are of order one. In this section we explore this region. 
The equations we employ are of course the equation of motion for "^rM modes (5.4) 
and the modified equation for the inflaton field (6.15). 

We follow the line of analysis performed in [30]. That is, in first step we ignore 
the AV term in (6.15), i.e. we study classical dynamics of the infiaton field when 
e is of order one and then study the dynamics of the \E'r,«m fields using this solution 
for (j) as background, and finally, we include AV in the equation of (p. The infiaton 
potential we start with (2.9) is a generic quartic potential and the analytic treatment 
of the equations is not possible for generic values of parameters Ae//, i^eff, ?77.^. The 
analysis for the chaotic case, i.e. when only Ae// (or when m^) is non-zero has been 
carried out in some detail in [15,30-32]. In our model, however, the case with only 
non-zero does not have the quartic coupling to the preheat field(s) and hence 
does not involve a preheating model. We therefore focus on the A0^/4 theory. 

For the A0'^/4 theory, the potential for the fields is 



V{<P, ^r,lm) = l^eff<P' + ^Ae//0^ E ^ " ^)Klm'^rM (6.19) 

r,lm 

1 1 
= TAe//0' + ^Ae//0'^|*C 



2 

Oml 



2 

m=l 

E E + E E 1*/^' 

/=0 m=l /=1 m=l 



where in the second line we have decomposed '^r,im into zero, a and /3 modes for which 
CO respectively takes values —1, —{I + 1) and I. We remark that the potential (6.19) 
is an approximation to the potential term we start with (2.1) (for the k — 0, — 
case) to order 

As shown in [30], the equation of motion toward the end of infiation is an 
(anharmonic) oscillator around = 0, the amplitude of the solution is decreasing 
as t^^/"^. In this regime the field is effectively living in a radiation dominated 
background with 1/H = 2t. It appears that the equation for takes a simple form 
in the conformal time, with Jacobi (elliptic) cosine function as its solution [30] . The 
equation for the modes, too, can be solved explicitly in our case. In fact in the 
notations of [30], the equation for all three zero, a and /3 modes is of the form 
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of g'^/X = n{n + l)/2 (with n = 1,1,1—1 respectively for zero, a and f3 modes) 
for which most of the calculations can be performed analytically. (Note that g'^/X 
of [30] coincides with the parameter u (6.8) in our model.) As discussed in [30] 
for these specific values of g^/X we have the significant property that there is an 
enhancement in the parametric resonance leading to considerable creation of zero, a 
and P modes. As discussed in [30] one can distinguish two even and odd n cases. For 
the odd n (i.e. for our zero modes, odd I a-mode, and even I /3-mode) the particle 
creation is peaked around zero momentum k modes. For the even n modes, however, 

the particle creations is peaked around momenta /c^ = \^inf^\[^' where e is the 
computed for the beginning of the slow-roll inflation and is the Hubble during 
inflation. For low n, the Floquet index /i^ oc Inn^ (n^ is the number density of 
the produced particles at momentum k) is around 0.15 for odd n and around 0.5 
for even n. Therefore, among the low n modes the main contribution to preheating 
is coming from odd n [30]. As discussed in [30] the bigger k is, the more energy 
can be transferred from the inflationary sector to the ^ sector and a more efficient 
preheat mechanism. This means that a and (5 modes with large I, I oi order N, 
make the biggest contribution, this is despite the fact that the zero modes have a 
larger degeneracy (of order iV^) compared to the degeneracy of order N for the large 
/ modes. All in all, due to the existence of the large / modes, and for large N in 
our model we expect to have a very efficient preheating model. The computations 
for the modes with large g'^/ X has been carried out in [30] and the only point which 
is different in our case is that their result should be multiplied with the degeneracy 
factor 21 + 1. 

As we argued we have an efficient preheat mechanism in our model. As a very 
crude estimate of the preheat temperature in our M-flation setup we may hence use 
an instant efficient preheating that all the energy of the inflaton field has gone to 
effectively massless zero modes by the end of inflation, leading to 

N^T^ ~ 3H^M^ , (6.20) 

where N"^ estimates the number of species and T is the preheat temperature. As we 
see, this is as if we have effectively an instant preheating model in which the maximum 
temperature achieved is lowered by 1/ -\/]V. As a rough estimate taking H saturating 
its current bound H ~ 10~^Mp and N ~ 10^ then the preheat temperature becomes 
of order T ~ 10^^ Gev. Reducing the preheat temperature to below GUT scale 
is in principle a positive feature, as it removes the problem with overproduction of 
gravitinos. 

7. Motivation from String Theory 

Here we argue that our M-flation setup presented in section 2 with non-commutative 
matrices and potential in the form of (2.2) is strongly motivated from string theory. 
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In the context of string theory, the world-volume theory of coincident p-branes 
is described by a (supersymmetric) U {N) gauge theory. In this system, the transverse 
positions of the branes, $/, / = p + 1, 9, which from the world-volume theory are 
scalars in the adjoint representation of U (N), are hence NxN matrices. For the case 
of our interest, p = 3, there are six such scalars. The DBI action for the system of N 
coincident D3-branes in the background RR six form flux (sourced by a distribution 
of D5-branes) is given by (e.g. see [33]) 

Here Is is the string scale and gs is the perturbative string coupling. The operator 
STr on a product of matrices is the trace of their symmetrized product. The induced 

metric on branes, Qah, is given by Qab = GMNdaX^d^X^ , where indicates the 
ten-dimensional positions of the branes and Gmn is the ten dimensional background 
metric. Here the indices /, J = 4, 5, ■ ■ ■ ,9 represent the coordinates perpendicular to 
the branes world- volume, the indices a, 6 = 0, 1, 2, 3 represent the brane world- volume 
coordinates and the capital letters M, = 0, 1, ■ ■ ■ ,9 indicate the ten-dimensional 
coordinates. The matrix is due to non-commutativity properties of the system 
given by 

g'' = 5" + ^[^',^'], (7.2) 

and Cfjo,,3 is a rank-6 antisymmetric Ramond-Ramond (RR) field which has two 
legs along the direction transverse to the D3-brane. 

We consider the ten-dimensional IIB supergravity background 

3 8 

ds"^ = -2dx+dx- - ^{x'f{dx^f + ^ dxidxi 

i=\ 1=1 (7.3) 

„ 2/t ^ 

where i,j indices, which are ranging over 1, 2, 3, parameterize three out of six trans- 
verse directions to D3-brane and x^ include three spatial directions along the brane 
and five of the transverse directions to D3-branes. With = Aglfi^ the above 
background, with constant dilaton, is a solution to supergravity equations of mo- 
tion. This background is very similar to the background of [34] (see also Appendix 
D of [35] for a discussion on the Matrix model on the above background). 

If we turn-off fluctuations along the directions transverse to the branes and the 
X* directions (this may be done if we compactify these three directions on a of 
very small radius), fix the light-cone gauge on the D3- branes, expand the action and 
keep up to order four in X^ , we obtain 
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(7.4) 



If we redefine 

X, = v/(27r)3(?, lj<^, (7.5) 

and upon addition of the four-dimensional Einstein gravity, the above action takes 
the form of (2.1) with the potential (2.2) once we identify the parameters as 

A = 2'Kgs , k = ^^== , ^w? . (7.6) 

Although from the brane theory viewpoint we need to choose rh^ and k such 
that (7.3) is a solution to supergravity, namely Xw? = 4^:^/9, since we presented the 
brane theory by the way of motivation, at the level of M-flation action we may relax 
this condition and take A, k and as independent parameters. It is also worth 
noting that Am^ = 4«;^/9 corresponds to the "symmetry breaking" inflation potential 
(3.1). Furthermore, the minimum (p — /i ior the potential (3.1) corresponds to the 
supersymmetric background where N D3-branes blow-up into a giant D5-brane. 

It is worth noting that in the brane theory setting the U (N) symmetry appears as 
a gauge symmetry, while in our model we took it to be a global symmetry. Promoting 
U{N) to a gauge symmetry does not change our analysis of the SU (2) sector and 
the corresponding inflationary dynamics. Due to the gauge symmetry, however, not 
all the modes are physical. Among them the zero modes can be removed by the 
gauge transformations and hence in the theory where U (N) is gauged only we deal 
with a and /3 modes. Thus in the gauge symmetry case, from the first A^^ -|- 1 iso- 
curvature modes in the mass tower, only one {I — a— mode) remains, which has an 
amplitude of order few percent of curvature spectrum. However this will not change 
the analysis of section 6.1, as the main contribution were coming from a or (3 modes 
of / in the window (6.10) which are also present in the gauged theory. The analysis 
of preheating mechanism of section 6 will remain valid because again zero modes do 
not have the main contribution. 

As a result of motion of D3-branes in the background C^^^ flux, two of the 
directions transverse to D3-branes blow-up into (fuzzy) two sphere, which in the large 
N limit behaves as a D5- brane with world- volume x S'^ [33]. In this geometric 
picture our inflaton field is nothing but the radius of this two-sphere. In this sense 
the effective infiaton field in our M-fiation scenario is closely related to the inflaton 
in the "giant inflaton" model of [36] (see also [37]). 

8. Discussion 

In this work we have presented a new inflationary scenario, the M-flation, in which 
inflation is driven by matrix valued scalar flelds. M-flation, hence, falls into the 
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general class of multi-field infiation models and shares positive features of N-fiation. 
Specifically, we used M-fiation to obtain a super-Planckian (large field) field variation 
during infiation. This leads to a considerable amount of gravity waves which can be 
detected in future gravity wave observations such as PLANCK [16,17], CMPOL [18] 
and QUIET [19]. Moreover, as we discussed, due to the scalings with powers of N, 
the dimension of the matrices, M-flation bears a solution to the fine-tuning problem 
of the couphng in the A0^/4 chaotic inflation. 

Here we focused on a special class of M-fiation scenarios with potentials of the 
form (2.2) where the potential is quadratic in powers of and their commutators. 
Within our three parameter family of the potentials there are interesting special 
models of infiation which were analyzed in section 3. One may, however, start with 
other forms for the potential. Specifically, if one starts from the string theory real- 
ization of the scenario, then from action (7.1) one obtains higher power corrections 
in terms of $i and [$j,$j]. Furthermore, the kinetic energy may also have a non- 
trivial form and one may combine the idea of M-flation with a DBI non-standard 
kinetic energy [39] . It would be interesting to see the spectrum of the adiabatic and 
iso-curvature perturbations for this case of "DBI M-fiation" . 

One of the observable effects of multi-field infiation models is that besides the 
usual power spectrum of the adiabatic fluctuations P-ji, we also have a non-zero power 
spectrum for the iso-curvature perturbation, Ps^i^- We analyzed the ratio Ps^i,„/Pn 
for various inflationary models up to the end of inflation. Inside the Hubble radius 
this ratio is close to unity but once the mode leaves the Hubble radius this ratio 
decays quickly towards the end of inflation. As shown, our analytical estimates of 
Ps^i^/Pn are in good agreement with the numerical results. In order to relate this 
ratio to the observed values from CMB, however, we should also supplement our 
M-flation scenario with a reheating mechanism. 

As we discussed in section 5.3, the iso-curvature perturbations '^r,im do not pro- 
duce entropy perturbations nor couple to curvature perturbations. This is due to our 
initial conditions resulting in the fact that they are classically frozen during inflation. 
This in turn implies that they do not carry energy up to leading order in perturbation 
theory. However, they can contribute to entropy perturbations at second order in 
perturbation theory or through preheating mechanism. Via the same mechanisms, 
the iso-curvature perturbations ^r,im can produce non-Gaussinities which are under 
intense observational investigations. 

As was discussed in section 6 our model naturally contains a preheating sector 
(essentially the ^^./m-Hiodes) which due to the large number of these modes works 
very efficiently, taking away the energy stored in the inflaton field. In order to 
complete our model, we need to have a reheating model via which the energy of 
the ^f-modes is transferred into the Standard Model particles. This is postponed 
to future works. As a possibility for reheating mechanism in string theory setup, 
where $i represents the collective positions of N D3-branes, we may imagine that 
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the Standard Model of particle physics are confined to the branes in the forms of 
open strings gauge fields A^i^\ The question of reheating would be how to transfer 
energy from the $i fields, more precisely from the '^r,im modes, to the open string 



In this work we have restricted the analysis to a particular solution where '^r,im 
fields are absent in classical inflationary dynamics. The inflaton fleld (f) is the pro- 
jection of $i along the N x N irreducible representation of SU{2), the Jj matrices. 
In the fleld space of — ^r,im. this corresponds to a straight inflationary background. 

In general, one may consider an arbitrary initial condition where \&r./m fields are 
turned on. In the field space of </> — '^r,im this gives a complicated curved infiation- 
ary trajectory. One such possibility with a more controlled dynamics is to take Jj 
to form a reducible N x N representation of SU (2) which consists of n irreducible 
blocks. In this case the classical infiationary dynamics of our theory reduces to that 
of n decoupled scalar fields, each with generic quartic potential. In this case the iso- 
curvature perturbations would be non-adiabatic and a significant amount of entropy 
perturbations can be created during inflation. Similarly, one expects a considerable 
amount of non-Gaussianities to be produced in this case. It would be interesting to 
build an specific model of M-fiation where ^r,/m fields are turned on during infiation 
and calculate non-Gaussianities and entropy perturbations produced and compare 
them with the observational bounds. 

Acknowledgment 

We thank Y. Farzan, N. Khosravi and K. Nozari for useful discussions. We thank 
K. Turzynski for some computational assistance. We specially thank B. Bassett and 
R. Brandenberger for comments on the draft and for many useful insights. H. F. 
would like to thank KITPC, and M.M.Sh-J the Abdus-Salam IGTP, for hospitality 
during the final stage of this work. A. A. is supported by NSERC of Canada and 
MCTR 

A. Symmetry breaking inflation 

Here we study inflation from the symmetry breaking potential in some details. Sup- 
pose inflation starts when 0j > The total number of e-folds N^, is obtained by 



where 0/ is the endpoint of inflation and for later convenience we have deflned x — 



gauge fields . 




(A.l) 
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(t)f{<pf — n) and y — (t)i{<pi — n) . As usual, define the slow-roll parameters 

where ' denotes derivative with respect to 0. Inflations ends when e = 1, which is 
used to fix x and 0/ 

X = 4M| + Mp^l6M|. + 2M|.//2 (A.3) 

The scalar spectral index at 0^ is — 1 = 2r) — 6e at (pi which can be used to 
eliminate y 



12 + ^144 + 8(1-71^)^ 



y _ 

Ml l-n-jz 



(A.4) 



Plugging these values for x and y in (A.l), wc find an equation for fi/Mp. Solving 
this equation numerically with = 60 and Us = 0.96 from WMAP5 central value, 
one obtains fi/Mp ~ 26. This in turn yields 0j ^ AAMp and 0/ ~ 28Mp. 

The COBE normalization, can be used to fix the value of Ag// 

Using 5h — 2 X 10^^ and the above values for (pi and /x, one obtains Ag// — 10^^^. 
This corresponds to A?" ~ 10^ as in chaotic inflation case. 

It it is also instructive to look into gravity wave amplitudes, determined by 
quantity r, defined as the ratio of gravitational perturbation amplitude to scalar 
perturbation amplitude at (pi. 

8/-. ^ 16,^2^" 

r = -(l-n,) + -M|- 
= 4(1 - n7^) + 32— ^ . (A.6) 

y 

Using the above values for (pi and one obtains r < 0.2 which is consistent with the 
upper bound r < 0.22 from WMAP5. 

The analysis when inflation takes place in regions ii/2 < (pi < /i and < (pi < ii/2 
is similar to the above. 
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